Neutron Interactions as Seen by A Segmented 
Germanium Detector 



I. Abt, A. Caldwell, K. Kroninger^, J. Liu*, X. Liu, 
^ ! B. Majorovits 

O ■ 

p^ , Max- Planck- Institut fiir Physik, Miinchen, Germany 

> 

o 

'^ Abstract 

The GERmanium Detector Array, GERDA, is designed for the search for "neutri- 
_ noless double beta decay" (0z^2/3) with germanium detectors enriched in Ge. An 

qj • 18-fold segmented prototype detector for GERDA Phase II was exposed to an AmBe 

neutron source to improve the understanding of neutron induced backgrounds. Neu- 
^ I tron interactions with the germanium isotopes themselves and in the surrounding 

^ ' materials were studied. Segment information is used to identify neutron induced 

peaks in the recorded energy spectra. 

The Geant4 based simulation package MaGe is used to simulate the experiment. 

.•^ . Though many photon peaks from germanium isotopes excited by neutrons are cor- 

V/-^ I rectly described by Geant4, some physics processes were identified as being incor- 

CN ■ rectly treated or even missing. 
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1 Introduction 



The GERDA (GERmaniuin Detector Array) experiment [1], designed for the 
search for "neutrinoless double beta decay" (Oz/2/3), is currently under con- 
struction in Hall A of the INFN Gran Sasso National Laboratory (LNGS), 
Italy. Neutrons produced near the germanium detectors by penetrating cosmic- 
ray muons can induce background events. In addition, neutrons from {a,n) 
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reactions in the surrounding rock are also a potential source of background. 
The study of neutron interactions with germanium isotopes as well as the 
surrounding materials is thus of great interest. 

Segmented germanium detectors will be used in GERDA Phase II. It has been 
shown that segment information is very useful to identify photon induced 
background [2]. It is interesting to check if segment information can also help 
with the identification of neutron induced background. 

In order to study the issues mentioned above, a GERDA Phase II 18-fold 
segmented prototype detector [3] was exposed to an AmBe neutron source. 
Energy spectra were recorded for each segment and the core. The segment 
information was used to identify peaks induced by neutron interactions. 

The Geant4 [4,5] based simulation package, MaGe [6], has been co-developed 
by the GERDA and MAJORANA [7] collaborations. The simulation of neu- 
tron interactions was verified by comparisons to data. 



2 Experimental Setup and Data Sets 



The detector used is the first segmented GERDA prototype detector. The true 
coaxial cylindrical crystal has a height of 70 mm and a diameter of 75 mm with 
a 10 mm hole in the center. It is 18-fold segmented with a 6-fold segmentation 
in the azimuthal angle and a 3-fold segmentation in the height z. It was 
operated in a conventional test cryostat. Details of the detector and its cryostat 
can be found in [3]. 

A 1.1 GBq isotropic AmBe neutron source was used in the experiment. The 
energy spectrum of the neutrons emitted from the ^Be(Q;, nf^G* nuclear reac- 
tion extends to 12 MeV. High resolution measurements of the neutron energy 
spectra of this kind of neutron source are presented in [8,9]. The dependence 
of the emittance of photons from the de-excitation of ^^C* on the neutron 
energy is described in [9]. 

The neutron source was located in a cylindrical paraffin collimator. The schematic 
experimental setup (not to scale) is shown in Fig. 1. The center of the col- 
limator was vertically aligned to the center of the detector and the distance 
between source and detector center was about 1 m. 

The core and segment electrodes were connected to charge sensitive pre- 
amplifiers. Their output was digitized using 14-bit ADCs in an XIA Pixie-4 
data acquisition system [10] with a sampling rate of 75 MHz, and recorded 
separately when the core was triggered. Two different gain factors were chosen 
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Fig. 1. Schematic experimental setup (not to scale). 

for four different measurements. The data sets are listed in Tab. 1. A low gain 
factor was chosen so that the energy range up to ~ 11 MeV could be covered. 
A high gain factor was chosen for measurements up to ~ 3 MeV. 

Two measurements were performed with the AmBe source present. They are 
referred to as HGdat (High Gain data) and LGdat (Low Gain data). In or- 
der to determine the background from the laboratory environment two more 
measurements without the source were performed. They are referred to as 
HGbg (High Gain background) and LGbg (Low Gain background). The data 
samples with different gains were combined for the study below 3 MeV. 

Table 1 

Data sets recorded with and without source. 





With AmBe Source 


Without AmBe Source 


DAQ Gain 


High Gain 


Low Gain 


High Gain 


Low Gain 


Emax [MeV] 


~3.5 


- 11 


~ 3.5 


~ 11 


No. of Events 


7.1 M 


4.7 M 


1.5 M 


4.7 M 


Name 


HGdat 


LGdat 


HGbg 


LGbg 



3 Core Spectra 



The total energy deposited in the germanium crystal was read out from the 
core electrode of the detector. Fig. 2 shows the core energy spectra for the 
data and background in the range of [0.08, 3] MeV. The thick line indicates 
the sum of HGdat and LGdat. The fine line represents the normalized sum of 
HGbg and LGbg. The trigger thresholds were set such that the spectra above 
100 keV were not affected. 



Eight photon peaks from the background were fitted with a Gaussian function 
plus a first order polynomial to normalize the background to the data. They 
are associated with the decays of ^^^Pb (352 keV), ^i^Bi (609 keV, 1120 keV, 



1764 keV, 2448 keV), ^28 Ac (911 keV), ^°K (1461 keV), and ^o^Tl (2615 keV). 
The numbers of events in the peaks determined by the fits were used to calcu- 
late the data to background ratios. The average of these ratios, 1.279 ± 0.003, 
was then used to scale the background spectrum. 

Fig. 3 shows the spectra in the range of [3, 10.2] MeV. In this energy range 
the background is small, as there are hardly any natural radioactive elements 
producing photons with such high energies. 



4 Neutron Interactions as Seen by the Core 



The main interaction mechanisms of neutrons with energies less than 12 MeV 
are thermal capture, inelastic and elastic scattering. Elastic scattering does not 
induce peaks that can be identified, because there is no photon emitted and 
the recoil energy distribution is too flat. Not only the production mechanism 
of the excited nucleus is important for the identification of a peak. The de- 
excitation mechanism has also to be taken into account. In most cases the 
nucleus de-excites instantaneously with the emission of one or more photons. 
However, it can also undergo internal conversion, in which case an electron 
from a lower shell is emitted instead of a photon. The excited nucleus can also 
be meta-stable and not de-excite instantaneously. 

Table 2 lists the processes identified in the core energy spectrum. If inelastic 
scattering happens inside the germanium crystal, the nuclear recoil energy is 
recorded as well as the energies from some of the prompt photons. In case 
of instantaneous de-excitation they are summed up: EineiasUc = E^ + Brecon- 
This causes an asymmetric peak with a long recoil tail on the high energy 
side. Internal conversions only create identifiable peaks, if they occur inside 
the crystal, otherwise the emitted electrons do not reach the detector. 

Table 2 

Type of neutron processes identified in the core energy spectrum. 



Production 


De-excitation 


Symbolic Notation 


Short Form 




instantaneous 


n+^Z^(^+i)Z+7 


-4z(n,7) 


thermal capture 


meta-stable 


(A+l)m2^(^+l)Z+7 


-4z(n,7™) 


inelastic 


instantaneous 


n+^Z^^Z+n' + 7 


^Z(n,n'7) 


scattering 


internal conversion 


n+^Ge^^Ge+ + n' + e" 


^Ge(n,n'e) 



Table 3, 4 list all the peaks observed in the core energy spectrum due to 
neutron interactions within the germanium crystal as well as within the sur- 
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Fig. 2. Core energy spectra with and without source. The normahzation procedure 
is described in the text. The energy range is [0.08, 3] MeV. Peaks induced by the 
AmBe source are indicated with their energies. 
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Fig. 3. Core energy spectra with and without source. The normahzation procedure 
is described in the text. The energy range is [3, 10.2] MeV. Peaks induced by the 
AmBe source are indicated with their energies. 

rounding materials: H, C, CI in the paraffin collimator, iron barrel of the 
collimator, Al, Ce in the aluminum vacuum can, and copper container of the 
pre-amplifiers. 

Pure photon peaks were fitted with a Gaussian function plus a first order 
polynomial to get the mean energies, FWHMs and the numbers of events in 
the peaks. The 596 keV peak from ^^Ge(n, n'7) does not have a Gaussian 
distribution. The treatment of this peak will be described in section 5. The 
662 keV peak associated with ^^'^Ce has a significant background contribution 
from ^^'^Cs. This was subtracted. The 692 keV peak from ^^Ge(n, n'e) is hard 
to fit because it is asymmetric and broad. It is also contaminated by other 
peaks nearby. The number of events in this peak was estimated by integration. 

The 4.4 MeV peak is due to photons from the de-excitation of ^^C* created in 
the AmBe source by the interaction, ^Be(Q;,n)^^C*. It is Doppler broadened 
because of the movement of the ^^C* nuclei. The width of this peak listed in 
Tab. 4 was determined by the fit. Above 5 MeV the widths of most of the 
peaks had to be fixed in the fitting procedure due to low statistics. The peaks 
that are not identified are marked with a question mark. 



5 Neutron Interactions as Seen by the Segments 



The energies deposited in each of the segments are read out separately. This 
provides more information about the interactions inside the germanium crystal 
than can be extracted from the core signal alone. 

For example, a photon with an energy of the order of one MeV has a mean 
free path of several centimeters in the germanium crystal. It most probably 
deposits energy in several different segments because of multiple Compton 
scattering. The result is a multi-segment event, in short MSE. In contrast, 
if there is only one segment with an energy deposition, it is called a single- 



Table 3 

Peaks observed in the core energy spectrum (see Fig. 2) due to neutron interactions. 



Fitted Energy 


Fitted FWHM 


Interaction Type 


No. of Events 


[keV] 


[keV] 






139.4 


1.6 ±0.2 


^4Ge(n,7™) 


3377 ± 520 


197.9 


1.9 ±0.2 


™Ge(n,7™) 


3306 ± 503 


499.8 


1.9 ±0.7 


70Ge(n,7) 


503 ± 186 


595.7 * 


- 


74Ge(n,n'7) 


(18.4 ±2.5) X 103 


662.0 t 


1.9 ±0.1 


i40Ce(n,7) 


2802 ± 188 


685.6 


1.4 ±0.2 


? 


628 ±111 


692 t 


- 


72Ge(n,n'e) 


~7000 


708.5 


2.4 ±0.5 


35Cl(n,7),36cl^36Ar 


782 ± 197 


721.9 


1.9 ±0.2 


? 


3502 ± 148 


843.4 


2.4 ±0.5 


27Al(n, n'7) 


1558 ± 202 


846.6 


2.4 ±0.2 


^^Fe{n, n'-f) 


2802 ± 196 


867.8 


1.9 ±0.5 


^3Ge(n,7) 


425 ± 129 


962.2 


2.4 ±0.2 


63Cu(n,n'7) 


1041 ± 129 


1014.3 


2.4 ±0.2 


27Al(n, n'7) 


1958 ± 123 


1164.1 


2.6 ±0.5 


35ci(n,7) 


646 ± 140 


1200.8 


2.8 ±0.2 


DEP of 2223 • 


2318 ± 122 


1326.9 


2.4 ±0.2 


63Cu(n,n'7) 


711 ±91 


1711.8 


3.8 ±0.1 


SEP of 2223 * 


5555 ± 133 


1778.9 


2.6 ±0.2 


27Al(n,7),28A1^28si 


469 ± 73 


2223.2 


3.8 ±0.1 


iH(n,7) 


79349 ± 300 



* The fitting of the 596 keV peak is described in a later section. 

f The background contribution to the 662 keV peak was subtracted. 

I The number of events in the 692 keV peak was determined by integration. 

•k SEP, DEP stand for Single Escape Peak and Double Escape Peak, respectively. 

? The unidentified peaks are marked with a question mark. 



segment event, in short SSE. The power of discrimination of MSE and SSE 
induced by photon using segmented germanium detectors has been shown in 

[2]. 



Table 4 

Peaks observed in the core energy spectrum (see Fig. 3) due to neutron interactions. 

Fitted Energy [keV] Fitted FWHM [keV] Interaction Type No. of Events 

3427 

3931 

4441 

4946 

6113 

6904 

7126 

7416 

7633 

7793 

7918 

f The widths were fixed during the fit. 

• SEP, DEP stand for Single Escape Peak and Double Escape Peak, respectively. 

? The unidentified peaks are marked with a question mark. 



85 ±7 


DEP of 4441 * 


2354 ± 263 


87 ±5 


SEP of 4441 • 


5873 ± 368 


92 ±2 


9Be(a,n)i2c* 


14672 ± 297 


4.9 ±1.4 


i2C(n,7) 


68 ±15 


7t 


35Cl(n,7) 


75 ±12 


7t 


SEP of 7416 • 


60 ±10 


7t 


? 


38 ±9 


7t 


35Cl(n,7) 


70 ±10 


7t 


56Fe(n,7) 


18 ±10 


7.1 ±2.1 


35Cl(n,7) 


21 ±8 


6.8 ±1.4 


63Cu(n,7) 


29 ±8 



5.1 Neutron Inelastic Scattering 



Compared to photon induced events, a special characteristic of neutron in- 
elastic scattering in a germanium crystal is that not only the photon energy, 
but also the recoil energy, is recorded. The inelastic scattering peak in the 
core energy spectrum has a high energy recoil tail and, hence, is much less 
significant than a pure photon peak with the same number of events, as seen 
in Fig. 2. However, it is possible to partially separate out the recoil energy 
distribution using information from the individual segments. The disentangled 
photon peak is much more significant than the original peak. This is due to 
the way a segmented detector can provide information about event topologies. 

Figure 4 shows the three types of events contributing to the inelastic scattering 
peak in the core spectrum. In all cases the scattered neutron escapes: 

(1) The nuclear recoil energy and the prompt photon energy are deposited 
in the same segment; 

(2) The nuclear recoil energy is deposited within one segment, the prompt 
photon deposits its energy in several other segments; 

(3) The nuclear recoil energy is deposited within one segment while the 
prompt photon deposits its total energy within another segment. 
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Fig. 4. Three topologies of neutron inelastic scattering inside a germanium crystal. 

In the first case, only one segment has a signal. The energies recorded by 
the core and the segment are the same, i.e. Ecore = Egeg = E^ + Erecou- 
Segmentation cannot help to disentangle the two energies. In the second case, 
the recoil energy can be observed in one segment. As the photon energy is 
shared between several segments, there is no peaked distribution in any single 
segment. This would partially be recovered by segment energy summation. In 
the third case, the recoil energy is observed in one segment, while the photon 
is observed in another segment. To disentangle the photon peak from the recoil 
energy distribution, energy spectra of all the segments are summed to get a 
spectrum of the energy deposited in any segment. In this spectrum, the type 1 
events produce the same distribution as in the core spectrum and type 2 events 
form a flat distribution. The type 3 events, however, create a sharp photon 
peak at the original energy and an enhancement in the low energy region, the 
recoil energy distribution. 

Figure 5 shows the any segment spectrum (fine line) together with the core 
spectrum (thick line) in the relevant energy ranges. Three peaks at 596 keV, 
834 keV and 1039 keV, associated with inelastic scattering, '^''Ge(n, n'7), '^^Ge(n, n'7) 
and ™Ge{n, n'7), respectively, are clearly visible in the any segment spectrum. 
The latter two are washed out in the core spectrum. 

It is possible to extract the number of each type of events in the 596 keV 
peak. The core spectrum was used to determine the total number of events, 
Nfotai- An exponential function was fitted to the shoulder of the peak asso- 
ciated with the nuclear recoil energy distribution. A Gaussian function was 
fitted to the 609 keV background photon peak on the shoulder simultaneously. 
The background below the recoil structure was obtained from interpolating 
the spectrum below and above the shoulder. The total number of events equals 
to the difference between the fitted exponential and the background. The er- 
ror was estimated by assuming different levels and shapes of the background. 
The fit is shown in the inset of Fig. 5. The dashed line represents the back- 
ground, the solid line shows the exponential plus the Gaussian function from 
the background peak. 
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Fig. 5. The any segment spectrum (fine line) and the core spectrum (thick line) in 
the relevant energy ranges. The inset shows a close-up of the 596 keV peak in the 
core spectrum. The fitting of recoil energy distribution is described in the text. 

The number of type 3 events, Ntypes, was obtained by fitting a Gaussian func- 
tion plus a first order polynomial to the 596 keV peak in the any segment 
spectrum. The small shoulder caused by the contamination with type 1 events 
does not change the results of the fit significantly. 

The determination of the number of type 1 events requires the study of the 
peaks purely induced by photons. This provides the probability that a de- 
excitation photon deposits its energy in exactly one or in multiple segments. 
The relative strength of the peaks in the core and any segment spectrum, that 
is, 7^(E^) = N,oreiE^)/Ke7iE^) = (NssE + Nmse) /NssE, directly translates 
to the relative rate of the total number of inelastic scattering events to the sum 
of type 1 and type 3 events, that is, 7^(E;"^'-*-) = NMai/iNtyp,, + Ntype^)- 



A Gaussian function plus a first order polynomial were fitted to eleven of 
the most prominent background photon induced peaks in the core and any 
segment spectra, respectively. The numbers of events in the peaks from the 
fits were used to calculate the ratio, TZ{E^). The points with error bars in Fig. 6 
represent the ratios calculated at different energies. A second order polynomial 
was fitted to get an estimate of the ratio at any energy, TZ{E). The number of 
type 1 events can be calculated as Ntypei = Ntotai/T^iE'"'''^''''^"') - Ntypes- 



The number of type 2 events can then be calculated as A^^ 



type2 



N^ 



total ' 



Ntypei — 



Ntypei- The results concerning event topologies in the 596 keV peak are listed 
in the second row of Tab. 5. The percentage of single-segment events, that is. 



Ntypei, out of the total number of events is P = Ntypei/ Nt 
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Fig. 6. The "core to any segment ratio" as a function of the energy. 

Table 5 

Numbers of events in the 596 keV, 834 keV and 1039 keV peaks. 

E [keV] Ntypel Ntype2 NtypeZ Ntotal 



595.8 


(1 ± 1) X 10^ 


(10 ± 3) X 10''^ 


7285 ± 218 


(18.4 ±2.5) X 10^ 


834.0 


[0, 380] 


[4100, 4700] 


2592 ± 186 


[6700, 7700] 


1039.2 


[0, 240] 


[2700, 3100] 


1429 ± 182 


[4100, 4800] 



The numbers of type 3 events in the 834 keV and 1039 keV peaks were ob- 
tained by fitting the any segment spectrum. Since there is no peak at these 
two energies in the core spectrum, it is impossible to get Ntotai from a fit. 
However, since the percentage V = Nfypei/ Ntotai decreases with energy, V at 
834 keV and 1039 keV should be less than ^(596 keV). Taking into account 
the relation, Nfypei = N total /T^ — Ntypes, the ranges of numbers of events in 
different topologies in the 834 keV and 1039 keV peaks were calculated. They 
are listed in Tab. 5 as well. 

The following steps were used to disentangle the recoil energy Erecou spectrum 
of inelastic scattering with a prompt photon of energy E^: 

(1) Exactly two segments having an energy deposition greater than 10 keV 
were required. 

(2) If one segment had an energy deposition in the range [Ej — 3a, E^ + 3a], 
where a was the detector energy resolution, the energy deposited in the 
other segment was used. 

The following steps produced the background to the recoil spectrum: 

(1) Exactly two segments having an energy deposition greater than 10 keV 
were required. 

(2) Two energy side-bands, [Ey — 6a, Ey — 3a] and [Ey + 3a, E^ + 6a] were 
defined. 

(3) If one segment had an energy deposition in the side-bands, the energy 
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deposited in the other segment was used. 

Fig. 7 shows the disentangled recoil spectra related to the 596 keV, 834 keV 
and 1039 keV photon peaks. The histograms start at 10 keV. The spectrum 
is dominated by the electronic noise below. The recoil spectra extending to 
~ 100 keV are clearly visible above the background. 
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Fig. 7. Recoil energy spectra corresponding to the inelastic neutron scattering with 
prompt photons of energies of 596 keV, 834 keV and 1039 keV. 



5.2 Internal Conversion 



If the excited state of a nucleus has the same spin as the ground state, internal 
conversion [11,12] is the predominant mode of the de-excitation. Since the 
mean free path of an electron emitted from internal conversion is about 1 mm 
in germanium, the energy of the electron and the recoil of the nucleus are 
deposited in the same segment. The core and the any segment spectra are 
the same. This is demonstrated in Fig. 5. The 692 keV peak from internal 
conversion, ^^Ge(n, n'e), is neither changed nor suppressed in the any segment 
spectrum. 



5.3 Double Escape Peaks 



The double escape peaks are enhanced in the any segment spectrum, because 
many events from the single escape and full energy peaks in the core spec- 
trum move to this peak. Two enhanced double escape peaks at 1200 keV and 
1592 keV are clearly visible in Fig. 5. They originate from the 2223 keV peak 
of iH(n, 7) and the 2614 keV peak of ^o^Tl. 
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6 Verification of Simulations 



MaGe, a C++ simulation package developed by the Monte Carlo groups of 
the Majorana and Gerda collaborations, was used to simulate the experiment. 
It is based on Geant4 [4,5]. The version Geant4 8.2 with patch-01 was used. 



6. 1 Generator, Geometry and Process 



Figure 5 in [8] shows the measured neutron spectrum emitted from an AmBe 
source. It was normalized to a probability density function and used in the 
neutron generator to assign energies to the outgoing neutrons. The genera- 
tor also produced 4.4 MeV photons from the ^^C* de-excitation inside the 
AmBe source. The Doppler broadening of the 4.4 MeV peak was simulated by 
Gaussian smearing with the observed widths taken from Tab. 4. 

The geometry of the experiment was implemented according to technical draw- 
ings. Approximations in the order of several centimeters had to be made re- 
garding 

• the shape and size of the AmBe source and how it is held inside the paraffin 
collimator, 

• the exact relative position between the crystal and the paraffin collimator, 

• the exact geometry of the components inside the cryostat. 

Geant4 provides high precision models for the simulation of interactions of 
neutrons with energy below 20 MeV [4,5]. The models depend on the "evalu- 
ated neutron data library" (G4NDL) for cross sections, angular distributions 
and final state information. The version G4NDL3.10 was used. 



6.2 Gore Spectrum 



Fig. 8 compares the core spectra of the measurement and the simulation in the 
range of [0.1, 3] MeV. The threshold effects below 100 keV were not taken into 
account in the simulation. The thick line is the experimental data. The fine line 
is the sum of the simulation and measured background. The background was 
normalized to data as described in section 3. The simulation was normalized 

to data according to the relation, Ndata = Nbackground + Nsignal = Nbackground + 

N simulation, where the Vs are the event numbers in the data, background and 
simulated spectra. Fig. 9 shows the same spectra in the range of [3, 10.2] MeV. 
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Fig. 8. Comparison of the neutron core energy spectra from 0.1 MeV to 3 MeV 
between data and simulation plus measured background. 
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Fig. 9. Comparison of the neutron core energy spectra from 3 MeV to 10.2 MeV 
between data and simulation plus measured background. 

6.3 Discrepancies between Data and Simulation 



The shapes of the continuous spectra from the simulation and data deviate 
due to the poor knowledge of the exact material and geometry of components 
between the source and the crystal. 

There is a known bug [14] in Geant4 concerning neutron inelastic scatterings. 
The secondary particles are not boosted back to the laboratory frame after 
the calculations in the center of mass frame are completed. This causes two 
problems: 

• The simulated recoil energies of the germanium isotopes are wrong. 

• The photon peaks from the interactions are not broadened. 

The first effect is demonstrated in the third inset of Fig. 8. The measured 
596 keV peak from '^^Ge(n, n'7) has a long tail on the high energy side due to 
the nuclear recoil, while the simulated peak misses this feature. 

The second effect is demonstrated in the inset of Fig. 9. The simulation gen- 
erates a broad and a narrow peak, both at 4.4 MeV. The broad peak is due 
to the de-excitation of ^^C* created in the source. The generator was adjusted 
according to the data. The narrow one is due to neutron inelastic scattering 
on carbon atoms in the paraffin collimator, ^^C(n, n'7). In reality, the carbon 
atom can gain a velocity of up to 0.02c causing a Doppler broadening of the 
order of 50 keV-100 keV. This is comparable to the broadening in the ^^C* 
de-excitation peak, and can, thus, not be resolved in the measured spectrum. 

The mean value from the Gaussian ffi to the measured 2223 keV photon peak 
is (2223.24 ± 0.01) keV. The simulated peak centers at (2224.61 ± 0.01) keV. 
This shifted value comes from the evaluated neutron data library. This problem 
has been reported to the Geant4 Problem Tracking System [15]. It was fixed 
for our studies by changing the value in the database to the measured energy. 
The result is shown in Fig. 10. 
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Fig. 10. The 2223 keV photon peak from H(n, 7) in data and simulation. The sim- 
ulated peak is shifted at 2224.6 keV before the modification described in the text. 

The 139 keV and 196 keV photon peaks from the meta-stable states of ^^Ge 
and ^^Ge produced by neutron captures are missing in the simulated neutron 
spectrum, see the first two insets of Fig. 8. This problem has been reported 
to the Geant4 Problem Tracking System [16]. 

The 692 keV peak from internal conversion, ^^Ge(n, n'e), is also missing in the 
simulation, see Fig. 8. It also has been reported [17]. 



7 Conclusion and Outlook 



An 18-fold segmented germanium detector was exposed to an AmBe neutron 
source and spectra were taken. A number of peaks from neutron interactions 
on germanium isotopes as well as the surrounding materials were identified. 
The segment information proved to be very helpful in identifying these peaks. 
Inelastic neutron scattering produces many events with energy depositions 
in more than one segment. Hence, the improved understanding of neutron 
induced interactions can also help to reduce the related background in the 
Oz/2/3 decay experiment, GERDA. 

The Geant4 based simulation package, MaGe, was used to simulate the ex- 
periment. Several discrepancies between data and MC were found. Further 
verification and improvement of the related Geant4 codes are needed. 

The experiment was not shielded from photons originating in the source. This 
resulted in a low signal (neutron induced peaks) to background (Compton 
shoulders of high energy photons) ratio. A new experiment including photon 
shields will be performed to suppress the photon induced background and 
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further investigate the neutron interactions on germanium isotopes. 
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